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Abstract 

In this paper we, to our knowledge, for the first time obtain the rate equations for Zeeman co- 
herences in the broad line approximation and steady-state balance equations directly from optical 
Bloch equations without the use of the perturbation theory. The broad line approximation allows 
us to use the adiabatic elimination procedure in order to eliminate the optical coherences from the 
optical Bloch equations, but the steady-state condition allows us to derive the balance equations 
straightforward. We compare our approach with the perturbation theory approach as given pre- 
viously and show that our approach is more flexible in analyzing various experiments. Meanwhile 
we also show the validity and limitations of the application of the rate equations in experiments 
with coherent atomic excitation, when either broad line approximation or steady-state conditions 
hold. Thus we have shown the basis for modeling the coherent atomic excitation experiments by 
using the relatively simple rate equations, provided that certain experimental conditions hold. 
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I. INTRODUCTION 



Coherent effects in laser radiation interaction with atoms and molecules play a major 
role in the physics and chemistry. Applications, such as electromagnetically induced trans- 
arency [l[ , laser cooling 0, 0, 0| , lasing without inversion , coherent population transfer 
different nonlinear magneto optical effects Q, new methods for magnetometry [§, 0], 
coherent control of chemical reactions and many other are widely used as a powerful 
research tools. Theoretical and experimental investigations of the coherent effects become 
increasingly important, as they open the way for more practical applications. Apart from 
some relatively simple cases, when direct solution of time dependant Scrodinger equation 
can be used |6[, usually, when one speaks about the modelling of experiments with atomic 
coherent excitation, he means the so-called "optical Bloch equations" (OBE), or Liouville 
equations 11,0, El for quantum density matrix p. These involve both optical and Zee- 



man coherences created in an ensemble of atoms. Zeeman coherences are quite stable and 
therefore it is relatively easy to employ them practically - for example, Zeeman coherences 
are a basic ingredient of sub-Doppler and sub- recoil laser cooling mechanisms jH, S 0] ■ Op- 
tical coherences, on the other hand, are very sensitive to a variety of factors - collisions, 
finite laser line-width, laser light fluctuations - both in phase and in amplitude, and many 
other. This means that, in describing a wide variety of atomic excitation experiments one 
can neglect the optical coherences. It leads to the well known rate equations for Zeeman 
coherences By saying "rate equations" we mean that they do not couple Zeeman 

coherences to optical coherences. 

Among the first ones to obtain the rate equations for Zeeman coherences by neglecting 
the optical coherences were C. Cohen- Tannoudji and J.P.Barrat in 1961 They used 

perturbation theory to obtain the rate equations in the so-called "broad- line approximation" 
(BLA) |l5||. These rate equations were obtained by considering the excitation with light from 
the spectral lamp and did not include the light induced transition effects into the analysis 
. They assumed intuitively that one can neglect the optical coherences in case of such an 
excitation. No mathematical arguments were provided, and the only justification of the 
used model was the good agreement between the theory and experiment. This lack of the 
rigorous mathematical argument was overcome later by C.Cohen- Tannoudji - with a slig htly 
different approach, through the use of the perturbation theory and assuming the BLA 1121. 
The BLA means, that the spectral line-width of the laser light used in excitation of atomic 
transition Au is very large compared to the natural line-width T of the atomic transition 

Acj>r, (i) 

and the spacing between laser modes 8ui is small compared to Y 

5u < T. (2) 

In this case different "Bennett holes" , burnt by the various modes in the Doppler profile, 
overlap, and the structure caused by different holes in the atomic response disappears. If, in 
addition, the modes cover all the velocity distribution, the atomic response does not depend 
on the velocity of translation motion of the atom and quantum density matrix p refers to 
internal variables only. In order to use the perturbation theory, the following condition must 
be satisfied: 

Acu>r,r p , (3) 
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where T p is connected to the time T p = — characterizing the evolution of density matrix 

p under the effect of the coupling with the light beam. The rate equations for Zeeman 
coherences, obtained by considering the conditions (jlj - (J3J), are often called the BLA 
equations. 

In the past rate equations for Zeeman coherences in BLA were very successfully used 
to analyze numerous nonlinear magneto optical effects. These include, for example, the 
interaction of molecules with multimode laser radiation - nonlinear Hanle effect, quantum 
beats, beat resonance, alignment to orientation conversion in a magnetic field etc., see, for 
example 3, 3 and references therein. This approach from the viewpoint of implementation 
in the form of computer routines seems to be substantially less demanding technically than 
the OBE. 

On the other hand, these rate equations for Zeeman coherences currently is not an often 
used approach to describe the laser radiation interaction with atomic gas. At first this seems 
obvious, as for the typical laser and atomic line-width the BLA conditions seem to be a very 
special case. 

However recently we have applied these equations for description of some linear and 
nonlinear magneto-optical effects in stationary interaction of alkali atoms with a broad 



band diode laser radiation |18l . I19L I20L |21[. In these cases the BLA clearly did not hold. 

Nevertheless, the agreement between simulation and experiment was good. The detailed 
analysis showed that for the use of BLA equations one does not always have to consider 
the BLA conditions - a rather striking result at a first sight. For example, in the case 
of a "steady-state" excitation there actually are no limitations for the use of the BLA 
equations except for the "steady-state" itself. The "steady-state" or stationary excitation 
means that the excitation light does not depend on time, which implies the same for the total 
density matrix p(t) - and this is the case in the large number of coherent atomic excitation 
experiments. 

What was the reason for such a good agreement between simulation and experiment in 
the above experiments with alkali atoms? After a detailed analysis it turned out, that the 
key factor was the fact, that the spectral line- width of radiation from the diode lasers was 
mainly determined by the phase fluctuations. The problem was, that the rigorous analysis of 
the limitations of the rate equations for Zeeman coherences in case of a noisy laser radiation 
seemed to be still lacking. On the other hand, there has been a large amount of work (see 
the overview in jiil . E^T ]) dealing with the OBE when the exciting radiation has a finite 
line-width arising from the fluctuations - both in phase and in amplitude. 

Thus in this paper we use the results obtained for the OBE and to our knowledge for the 
first time obtain the rate equations for Zeeman coherences dir ectly from the OBE. We also 
compare our approach with the perturbation theory approach |l2| and show the advantage 
of our approach. We analyze the limitations of usage of the rate equations for Zeeman 
coherences in conditions of noisy laser radiation - with an accent on analysis for nonlinear 
magneto-optical effects in atoms. In the limit of large angular momentum (molecular case) 
such an analysis, at least partially, was done previously in In this paper our goal is 
to fill in this gap for the atoms. The obtained results are in such a good agreement with 
experiment, that we feel, that the relatively simple rate equations' approach (comparing to 
the conventional OBE approach) is far too often undeservedly neglected, when discussing 
the modelling of nonlinear magneto-optical effects in atoms. 
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II. BROAD BAND RADIATION INTERACTION WITH ATOMS. 



A. Exciting light. 

In our analysis of usage of rate equation for Zeeman coherences we will describe the 
exciting light classically by a fluctuating electric field E(i) polarized along the unit vector e 

E(t) = e(t)e + e*(t)e*, (4) 



e (t) = \ £u \ e -<*(*)-i(«-ktov)t. (5) 

We account for a shift uJ — k^v in the laser frequency due to the Doppler effect - v is the 
velocity of translation motion of atoms and kjj is the wave vector of the exciting light, uj 
is the center frequency of the spectrum, |e^| is an amplitude of laser light field and $(£) 
is the fluctuating phase, which gives the spectrum of the radiation a finite bandwidth Auj. 
If the phase fluctuations are completely random, then the line-shape of the exciting light 
is Lorentzian. In the case of a laser this light corresponds to the single mode laser with 
randomly fluctuating phase. In the case of a spectral lamp this light corresponds to the 
lamp, where the dominant mechanism for the line- width broadening is determined by the 
collisions between the radiating atoms or molecules. Note, that for the single-mode laser 
the BLA condition (j2J is not fulfilled. 

The Rabi frequency is determined by 

n* = (6) 

where d is assumed to be the strongest atomic electric dipole moment for the transition 
(transitions) under consideration. 



B. Optical Bloch equations. 

We consider the dipole interaction of an atom with a laser field in presence of an external 
static magnetic field B. We assume that the atomic center of mass moves classically, which 
means, that the only effect of the dipole interaction of the atom with a laser field is an 
excitation of a classically moving atom at the internal transitions. In this case the internal 
atomic dynamics is described by the semiclassical atomic density matrix p, which paramet- 
rically depends on the classical coordinates of the atomic center of mass. We consider atoms 
with definite velocity v, illuminated by the exciting light ({! )) . (|5 ]L resonant with the g <-> e 
transition, in presence of an external static magnetic field B, which removes the degeneracy 
of the levels g and e, so that now we consider Zeeman sublevels g± and e^. In writing OBE, 
see for example [22] 

+ ihRp, (7) 

we consider only the relaxation R due to spontaneous emission. This means that we neglect 
other relaxation mechanisms, such as collisions, fly-through relaxation etc. This assumption 
means, that different velocity groups do not interact - the density of atoms is sufficiently 
low. For simplicity we also assume that the atomic transition forms a closed system - cycling 



op 



H,p 
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transition. In this case the spontaneous relaxation terms for a closed system for the density 
matrix elements p gigj , p giej , p ei9j , p etCj are: 



R P9r9j 


9r9oP^ e 3 








1 


R P 9i ej 






r 


RPe l9j 


~~ ~~^Pz%9ji 




= Peiej i 



(8) 



where T e g f gj describes the spontaneous relaxation from p eiej to p gigj and T describes the 
spontaneous relaxation from e — > g. For the closed system it is obvious that J2^7i9j = T- 

9i9j 

Hamiltonian H = H + V includes the unperturbed atomic Hamiltonian H , which depends 
on the internal atomic coordinates: Hq \ ^ n ) = E n \^ n ), and the dipole interaction operator 
V = — d • E(£), where d is the electric dipole operator. Writing OBE explicitly for the 
density matrix element pij, we get: 



dp, 



dt 



i 
h 



H , Pi. 



h 



d-E(t),pJ +Rp 



= -iuJijPij + Rpij + 



(9) 



h 



n k n k 



Ei — Ej 

where uj^ = — - denotes the Zeeman splitting of the levels i and j. By choosing 

it 

quantization axis to be parallel to the external static magnetic field B, all the dependence 
of the density matrix on the B field is included in the splitting term lo^. Thus we arrive to 
the following equations for the density matrix elements p gigj , p giej , Pe, gj , Pe iej '- 



dp, 



9i9j 



dt 



h 



(9i\d-e\e k ) p ek9j + ^e* |d • e*| e fc ) p ek9j - 



efe 



efe 



^2 (e fe |d • e| 9j ) p giCk - \e*Y^ (e k |d • e*| 9j ) p gi£k - iuj gigj p 9igj + ^Tl^'p. 



9i9j' ' 



(10) 



dp gi£j i ST / \A \ \ 

" T£2_^\9i\d- e\e k ) p eke 



dt h 



+ h £ *Yl (9i\ d - e *\e k ) p ekCj 



efe 



efe 



~ Ji £ Yl l d ' e l e i) P9i9k - l £ *Yl ( 9k l d ' e 1 e i) PtoQk - iUgiejPg iej ~ \pg^ , (11) 



9k 



h 



9 k 
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<J.t 



()f jfY ( ei l d ' e l g k) Paw + fr £ *Y ( ei l d ' e *l 9k) Pg. 

9k 9k 



e-k e.k 

= jfY ( &i l d ' e l 9k) P9kej + J^Y ( 6i l d ' ^1 9k) Pa^e- 
9k 9k 

~ jfY ^ 9k ' d ' G l 6 ^ Pe ^ ~ jfY ( gk l d ' ^1 e ^ P*iOH ~ '^■<,P<>, - Fp eiej , (13) 
9k 9k 

The matrix elements of the type (e*| d-e \gj) can be calculated using the standard angular 
momentum algebra 12 , 2^, Hf| . 

Now, in order to eliminate the fast oscillations with optical frequency 57, we make the 
following substitutions: 

P9i9j = P9%9] = P9i9ji 

P 9iej = P^e/ { ^ )t+m \ (14) 
n - n p -i(p-*3Y)t-i$(t) 

P&i^j P&i&j P&i&j' 



By using the rotating wave approximation [11| and neglecting terms with double optical 
frequency, we arrive to: 



dp % ^ ^ ^ 

~i}T = n \ ejj \ Y i d ■ e *i ek ) p°i*i ~ n Y ( ek i d ■ e i 9j) p^k- 

~ '-;/<;/.,/'•/ ./. + Yj'fnfuP'-'. ' S '•">! 



dp. 



9i<-'j 



Of h ^ fa l d • e*| e k ) p ekej - - \ejs\ Y (9k |d ■ e*| e,-) p 9l9fc - 

efc ' 9k 

_ _ r <9$(t) 

- *(W - k ^ v + ^j)Pmi - g Put* ~ i—fa-PfH*i ' ( 16 ) 

^If 1 = ^ N ^ (ei \d ■ e| </ fc ) p gs9 . - - |^| ^ (e k |d • e| p ei£k + 

9k e fe 

.,_ . , r .d$(t) 

+ i(w - k^v - u eigj )p eig . - - Pei9j + i—^-p ei9j , (17) 
dp % ^ ^ ^ 

-^p = ^ m i d ■ e i ^ - ^ i d ■ e *i e ^ 

- iu^p^e, - Tp eie (18) 
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C. Atoms in a fluctuating optical field. 



The equations ()15|) - (fTHjl are stochastic differential equations |27( with stochastic variable 
<9$(t) 

— - — . In an experiment, as a rule, we deal with quantities that are averaged over the time 
at 

intervals that are large in comparison with phase fluctuation time in the excitation light 
source, therefore we need to perform the statistical averaging of the above equations. In order 
to do that, we solve the equations lfT^j) and (|T7|) (with initial conditionp 9 . ej . (t ) = p ei9j (t ) = 0) 
and then take a formal statistical average over the fluctuating phases: 

= ^ kd E l d ■ e *l e k) (pZg 3 ) - kd E ( Ck l d ' e l to) (p^ k )- 



tuJ gigj 



(/w) • £n$ow> ( 19 ) 



t 

= 'ft kd E & l d ■ e *l e ^ J el-^-^^)-^') ( Pefcej (t')e- 4[ * W ^ (t,)1 ) dt'- 



t 

I 



- -ft kd E (9k Id ■ e*| e 3 ) J el-^-^+^-^t-t') ^(fle^®-*?*) *! , (20) 

9k to 

t 

fe) = { kd E ^ > d ■ e > 9k) I e^-^-«W-m-V (^(^») dt- 

9k tn 



to 
t 



- -ft kd E < e * l d ■ e l 9i) J eK-~ k — (^(f)^ 01 )*', (21) 

to 

Q /p \ ^ ^ % ^ 

— = ^ 1^1 E ( ei i d ■ e *i ^ (piw) - % kd E ( gk i d ■ e *i e ^ (^»)~ 

- izu^e, (p eie . ) - T (p eie . ) . (22) 

Now we employ the relation (JSJ) which allows us to use the decorrelation approximation 
[iii liol The decorrelation approximation means that we neglect the fluctuations of 
Pa ia:j {t) {cl = e,g) around their mean value (p^o,^)), and thus separate atom and field 
variables in JUJ) and (J2U): 

(pa.aAt> ±im) ~* {t ' )] ) = (P^W) (e^-^) , (23) 

where a = e,g. The decorrelation approximation in general is valid only for Wiener-Levy- 
type (see below) phase fluctuations [3(1 Hl[. In the case of a general stochastic field the 
decorrelation approximation can be used as a first approximation only for weak fields below 
saturation 29l 1301 . 
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In order to evaluate the correlation function ^e ±l t* < -*- ) ~ <I> (''^), we assume two sim ple models, 
which lead to similar results. The first one is the "phase jump" model (3^, I33I. l34| . which 
assumes the phase to remain constant, except sudden random "jumps", when it changes 
to a new constant value. This model is used to describe the fluctuations from the spectral 
lamp, when the dominant mechanism for the line- width broadening is determined by the 
collisions between the radiating atoms or molecules (221 . I32L Isl^ . The second model is the 
"phase diffusion" model 3(1 31 , |35|, |36| , which assumes the continuous random diffusion of 
the phase. This model is used to describe the fluctuations from the single-mode laser with 



fluctuating phase ji^, IH 36| 



1. "Phase jump" model. 

Our analysis of phase jumps in excitation radiation is based on the detailed analysis of 



the model in case of optical Bloch equations performed in [32J, |33|, |3J]. In this work we 
closely follow a n ap proach, that is analogous to the effect of instantaneous collisions on the 
density matrix 22]. The random jump process is Poissonian in nature - the probability for 
the phase to change N times during time t — if is 



1 ft -if 



N 



P «=MVJ (24) 

where T is the average time between successive phase jumps. Now we define (|e ±lA *) as the 
average phase change during one jump. If during the time t — if there has been only one 
phase jump, then (er ll ' l ^^~^ t = ^e ±lA *). Obviously, if during the time t — t' there has 

been N phase jumps, then ^e^*^ - *^)^ = (e ±lA *) iV , as every jump in average adds one 
more multiplier ^e ±iA<I> ). In order to get the final expression for ^e ±l '*^ - *^'^), we must 
average over every possible number iV = 4- oo of phase jumps during time t — if: 

OO OO 1 , 4-4-1 \ N 

JV=0 " A N=0^' ^ ^ ' 

= e-^( 1 -( e±!A *)). (25) 

At this point we make further simplifications. We consider the case, when there is no 
correlation of the phase values before and after the jump Then T is also the correlation 

2 

time of the phase (in average after the time T the phase "forgets" its past) T = . We 

Au 

also consider, that all the phase values occur with equal probability. Then (|e ±lA *) = 0, and 
fl23|) becomes 

^ e ± W )-$ (t ')]\ ) = e -^( t - t ')_ (26) 

This means that the spectral distribution of the exciting light is Lorentzian with FWHM 
Au. 



2. "Phase diffusion" model. 

Our analysis of the influence of phase diffusion of the excitation radiation on the in- 
teraction of laser radiation with atoms we base on the phase diffusion model analyzed in 



S 



|29l l30L l31l . l35l |36(. In phase diffusion model the field has a constant amplitude, but its 
phase is a fluctuating quantity, which obeys the Langevin equation 



dt 



where q(t) is a Gaussian random force with correlation function 

(,{tW)) = bPe-^ t - t \ 



(27) 



(28) 



(?(*)> = 0, (29) 
which means, that q(t) obeys the Langevin equation for Brownian motion 2^, 22 



d 



(30) 



where F(t) is a 5-correlated Gaussian force fulfilling the correlation 

(F(t)F(t')) = 2b{3 2 5(t-t'). (31) 

The meaning of the parameters b and (3 can be interpreted from the equations (JTTj) - 

(13 — is the correlation time of the phase time derivative s(t), but b gives the band- width 

of the field in the limit /3 — > oo. Explicit expressions for /3 and 6 in terms of fundamental 
laser constants are discussed by Haken in |37J, see also [38j. 

The spectrum of the exciting radiation, described by (}2Tj) - (}3Tf . is given by the Fourier 
transform of the correlation function 



,±<[*(t)-*(f)] 



cxp 



/5 l 



- 1 



(32) 



For f3 ^> b the spectrum is Lorentzian with FWHM Auj = 2b and having a cut-off at 
frequencies (3, but for f3 -C b the spectrum is Gaussian with FWHM a/8 ln(2)6/?. 

In the limit /3 — > oo the spectrum is pure Lorentzian with FWHM Auj = 2b and s(t) 
becomes 5-correlated 

(<;(t)<;(t')) = 2b6(t-t'), (33) 

but the phase $(t) obeys a Wiener-Levy stochastic process. As mentioned above, the 
Wiener-Levy process is the only one, for which the decorrelation approximation is mathemat- 
ically rigorous [1^,1^, 0]. It is easily understood, as in this process the relevant fluctuating 

quantity ^(t) is 5-correlated (correlation time — of q(t) tends to zero when f3 







oo 



and 



thus we can always separate the time-scales of evolution of (pcna^t 1 )) and ^e ± *'* < -*' )-<I> ^'^) 
in (j2HJ). Wiener-Levy stochastic process is a nonstationary Markov Gaussian process 
and is described by the Langevin equation for Brownian motion with negligible acceleration 



27 , 130L HH, which can be shown to be equivalent to the diffusion equation 22|| . For 
Wiener-Levy process the relation (J3~2*|) becomes: 



,±i[*(t)-*(f) 



exp 



[-b\t-t'\ 



i(t-f) 



(34) 
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where we have used the fact, that t ^ t'. 

The Lorentz profile is not a good description of the wings of any laser spectrum, thus this 
model is appropriate for rather small detunings. However, as shown in (35l . 36l |. for ft ^> b 



(or ft ^> Au) in the limit ft ^$>T, Qr, the line-shape of the exciting light is Lorentzian with 

Au; 

cut-off at frequencies ft. In this case the damping term — — is simply multiplied by the 

cut-off term, dependent on the detuning |35l |36| . This corresponds to a more realistic model 
of the laser spectrum. We also have to remember, that for the rotating wave approximation 

d$(t) 

uoq 3> ft, as ft is the correlation time of the phase time derivative <;(t) = — — — . 

D. The effective relaxation caused by the fluctuations of the exciting light. 

As can be seen, both approaches give similar results for time average phase fluctuation 
value - the effect of the phase fluctuations on the density matrix is simply to add the 
additional relaxation term, equal to the HWHM of the exciting light. Now we use (|2fi|) and 
to rewrite (|19|) - (|22p (for simplicity we further drop the averaging brackets): 



~oT = h ^ E & l d • e *l e k) p^si ~ i M Yl ( ek l d ' e l 9j) 

efc e fe 



efc 



P^+jyT^, (35) 



C i (:' j 



t 



to 
t 



l - \e*\ (9k Id • e*| e,) J eH(^ k ^-,)-(^^)]^^p^(t')^, (36) 

9k t 



t 

. I 

Pei9j = 

9k 



£ (e* Id • e| g k ) j eK^-^^-C^l^p^CO^- 



to 

t 

i 



h 



(en Id • e| g s ) J e K^ k — ^hi^)}^ Peiek (t')dt', (37) 

to 



dpe t 



dt H 



|e«| E ( 6i l d ■ e \9k) 



h 



\e<jj\^(g k \&- e*\ej) Pei9k - 



9k 



9 k 



%UJ g^ g . p g^ g . ^^/^ 6 j £ 



(38) 



10 



III. RATE EQUATIONS. 

For the sake of simplicity we further assume (with ujb characterizing the Zeeman 
splitting): 

Au;>u; s , (39) 

though this condition can be avoided at the expense of complication of the final rate equa- 
tions. (|39|) means, that we can write: 

[V AuA ._ , fV Auo\ ,_ , 



/T AuA , , , /r AuA 



-i (a - k^v - u ) . (41) 



At this point we go further and assume certain conditions, which allows us to simplify 
significantly the expressions for optical coherences (}36|) and (|37|) . These conditions are either 
BLA (P) - (jHJ) or the steady-state ( (fHOjl - see below) conditions. Under these circumstances 
the expressions for optical coherences ({BE} and (pT7j) become: 



Piwy = ~k (t i Auj\ ~ r r ( & l d " e *l efc ) P*w ~ Yl |d • e*| ei ) p 9 . 

ft l2+T)+H w - k-v - uV 4 



(42) 



(43) 

Now, by substituting (jj2j) and in and (f3~%J) we arrive to the final rate equations: 



dp 



9i9j 



Of 



: T P Y ' dl ' ^1 6k ^ ( &m ' dl ' e ' Si) Pe k e m ~ 

(y + «AE P j (9i l d i ' e *l e k) (e* |di • e| # m ) p 9mffj - 

e fe ,g m 

( y - «AE P J (# TO |di • e*| e k ) (e k |di ■ e| p 9l9m - 

e fc ,g m 



(44) 



<9t 



T p Y ( e * l d i ' e l 9k) (9m |di ■ e*| e,) p 9fc9m - 

9fc.3m 

^y - iAE^j ^ (ei |d 2 • e| g k ) (g k |di • e*| e m ) 



(45) 



fy + «AEpj ^ (e m |di • e\g k ) {g k |di • e*| e,,)p e . e 

^ ' 9fe,e m 
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where d± = — denotes the electric dipole moment unity vector, and thus matrix elements 
|d| 

(ej| d • e \gj) are written as: 

(dl d ■ e \g s ) = (e;| d x • e \g s ) {e\\d\\g) , (46) 

where (e||d||<7) is the so-called reduced dipole matrix element. Note, that for the steady-state 
situation we must consider condition (JoTJj) in the above equations. T p and AE p are defined 
as: 

r \f i 2 f- + —) 



h 2 



x |(e||%)r x ^ 2 2 -, (47) 

+ +(^-k^v-^ ) 2 

A^=i^ixi W ii g )i'x , r+ 4r k !rr° — f w 

" (2 + —J + (<*>- k^jv - ^ ) 2 

T p is the probability per unit time of an absorption or stimulated emission process, and 
AE P describes the light shifts |l4| produced by the light irradiation (dynamic Stark shift). 
For BLA conditions © - © flll|) and (@HI) become: 

' " ~ '"^ >< l< e IMIs>| 2 x —a ; ^ -, (49) 



ft2 /A 



+ (uJ - k^v - u; ) 

A£ p « x | (e „ %) | a x f^" (50) 
ra (— J + (w - k^v - u; ) 2 

Note also, that the phase fluctuations (described by the above models) reduce the satu- 

_ r 

ration on resonance (a; — k^v — u>o = 0) by the factor — — , and increase the saturation 

r + Au 
r ~\~ /\ to 

far-off resonance (u — k^v — u 3> T, Au) by the factor . 

When the density matrix for the excited state is calculated, one can obtain fluorescence 
intensity with specific polarization along the unit vector ei as 

1 (ei) = To^2 ^ l d ' e *l 9i) (9i |d • ei\ej) p^, (51) 
where J is a proportionality coefficient. 



IV. ANALYSIS AND CONCLUSIONS. 
A. Perturbation theory approach. 

The obtained rate equations for Zeeman coherences coincide with equations obtained 
earlier in perturbation theory approach in ^2|. In perturbation theory approach as small 
parameters are used ratios between rate constants involved in the problem (T p , T) and line- 
width of the excitation radiation Au. Here we would like to stress that, however the obtained 
equations coincide, the approach used in this study is different and allows us to examine in 
more detail the limits of usage of rate equations for Zeeman coherences to analyze specific 
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experiments. To compare both approaches, let us have a brief look in method used and 
conclusions obtained with perturbation theory. 

Let T p = — be the time characterizing the evolution of density matrix p under the effect 
1 p 

of the coupling with the light beam. In the following analysis it is assumed that the intensity 

is sufficiently low so that T v is much longer than the correlation time T = — — of the light 

Auj 

wave 

Auj > r„. (52) 
Now consider a time interval At such that 

T p , r > At > T, (53) 

where r = — . Since T p , r 3> At, one can conclude that p(t+ At) —p(t) is very small and can 

be calculated by perturbation theory. By using perturbation theory, it is shown in 2| , that 
the average variation of p, (p(t + At) — p(t)) (the average is taken over all possible values 
of the random function e(t) - see below) is linear in At and only depends on p(t) 

(p(t + At)-p(t)) Ap(t) 

At = "AT" (54) 

This means, that we can replace — with the time derivative ^[ - , provided that 

1 At dt ' 1 

we never use - to describe the changes of p(t) over time intervals that are shorter than 

correlation time T of the light wave, which drives the atoms. — ^— = ^[ - is called the 

"coarse grained" derivative [40]. 

In jl2j the exciting light is taken to be the superposition of parallel plane waves having 
all the same polarization e, but different amplitudes \e^\, frequencies and phases $ M 

E(t) = e(t)e + e*(t)e*, (55) 



e(t) = I £ mI e-**"-*^"-^*. (56) 

BLA relations ((TJ - (0) hold and the relative phases of the different modes are assumed 
to be completely random and thus obeying the correlation relation: 

(e-«*«-V>) = V. ( 57 ) 

The instantaneous electric field e(t) of the light wave thus may be considered as a sta- 
tionary random function, which obeys the correlation relation: 



(e(t)s*(t - T )) = J2 1^1 Ivl (e-^-V)) e -i{^-^)t e i{^-^){t-r) = J- |e M |- ( 



2 -i(w^-k— v)t 



(58) 

Applying perturbation theory, after some calculations with the consideration of (|58j) and 
'coarse grained" derivative, the rate equations are obtained |l2j, again considering condition 
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(|39|) for simplicity. The obtained rate equations are exactly the same as the above derived 
equations flUI , ()45|) . but with T p and AE p having slightly different form as defined in ()47|) 

- (foUj) . This mismatch is easily avoided, if instead of the exciting light model ()55|) - ()58|) we 
take the model (IH), (|SJ), (J2E1), (jHlj) - Then the rate equations are the same as equations (JHJ), 
(JUJ), with T p and A£ p defined as in (J4l3 |) and (f5n j) . 

Thus the perturbation theory approach is summarized as follows: we define " a priori" the 
BLA conditions (0) - Q and then use the perturbation theory to obtain the rate equations 

- and thus we are restricted to the BLA case. 

However, in our approach we obtain the "phase- averaged" OBE and then it is possible 
to choose between the BLA (JH) - @ or steady-state (JoTij) possibilities. 

Thus we arrive to the conclusion stated above that the approach discussed in this article 
allows to examine the limits of usage of rate equations for Zeeman coherences to a greater 
detail than the perturbation theory approach. Therefore it can be applied to analyze larger 
number of experimental situations. 



B. Velocity dependence. 



When we look at (|4*Tj) - (fHUj) . we see, that T p (induced transition rate) and AE P (dynamic 
Stark shift) are velocity dependent, and thus are also the equations (JSJ), (|43j) . This means, 
that in describing the observable signal we need to take into account all the velocity groups 
involved (note, that we have already assumed, that different velocity groups do not interact - 
the density of atoms is sufficiently low). In standard method one has to determine the signal 
dependence on velocity and then sum (integrate) over the velocities (of course, assuming that 
velocity distribution is known). However, usually the signal dependence on velocity cannot 
be found in analytical form, as can be seen from (}4*4*|) - (|5Uj) . Thus a large amount of 
calculation is necessary to determine this dependence - and still it is just an approximation. 

The situation is simplified only for a specific kind of experiments. For example, if we 
consider a case, when exciting line- width Auj is much larger than Doppler width Aujd of 
the atomic line (as it was originally assumed in the perturbation theory approach as given 
inQ) 



Alu »Awd, (59) 

then, as mentioned above, the atomic response does not depend on the velocity of translation 
motion of the atom and quantum density matrix p refers to internal variables only. In such 
a case we obtain the rate equations by simply putting k^v = in (J4*4|) . (jUJ), which is the 
same as to consider the atomic velocity group kjjv = only. Only one velocity group is also 
involved in experiments with cold atomic gases, atomic beams, etc. 

However, we have successfully used the rate equations for Zeeman coherences (JHJ), (|43|) 
in modeling of various experiments 3 3 HI 21 1. In these experiments ()59|) clearly did not 



hold, nevertheless in describing experimental signal from all velocity groups, we have used 
the calculated signal from just one velocity group k-^v = (note that AE p (k^v = 0) = 0). 
It is clear, that in this case for the experimental and simulation results to coincide, we cannot 
use the exact expressions (}4Tj) . (j4T?j) for T p (k^v = 0). Thus we must consider the "effective" 
induced transition rate T p ^ , which in general does not coincide with T p (k^v = 0). 

Using the signal from velocity group k^v = as the calculated signal is justified, if we 
know the relation between T p (k^v = 0) and T p ^ in advance. In reality this relation is 
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known only in some specific cases - for example, for the " steady-state" excitation with laser 
intensities below saturation - then we know, that T e Jf ~ T p (k^v = 0), as the signal from 
velocity group kjjv = is proportional to the signal from all velocity groups - see fl^ . 
However, in most cases establishing the relation between T p (k^v = 0) and Ty* is rather 
complicated, as it involves a large amount of calculations. 

Therefore in analysis of experiments we have used the following approach - the signal 
from velocity group k^v = is calculated and then the best fit to an experiment is found - 
thus experimentally finding the relation between T p (k^v = 0) and Tff*. In order to predict 
further results, we use the extrapolation and various other mathematical techniques. This 
method has proven to be successful in many cases. 



C. Steady-state excitation. 

As it was shown above generally, the usage of the rate equations for Zeeman coherences for 
description of time dependent behavior of atoms in laser and magnetic fields requires certain 
conditions regarding absorption rate connected with light intensity and spectral width of 
the laser line. At the same time very often in coherent atomic excitation experiments the 
"steady-state" or stationary excitation conditions are reached - the excitation light does not 
depend on time, which implies the same for the total density matrix p(t). For the steady- 
state to happen, the system has to go over a rather large number of cycles, which means, 
that the steady-state is reached only after some time, after which it remains in this constant 
state forever (unless, of course, the conditions imposed on the system are changed). This 
means, that mathematically we can obtain the steady-state, if we consider p(t = oo) at some 
certain time moment in eternity t = oo, when we can be sure, that the system has reached 
its steady-state - of course, if it can reach the steady-state at all. It is also obvious, that 
for steady-state the time derivative of the density matrix is zero and thus mathematically 
we can also obtain the steady-state by simply putting dp(t)/dt = for both optical and 
Zeeman coherences. 

Thus under the steady-state conditions we can express the optical coherences in terms of 
Zeeman coherences from the OBE straightforward, without any assumptions. In doing so 
we obtain the rate equations for Zeeman coherences p 9igj {t) and p eiej (t), which now form a 
set of linear equations, because of the steady-state condition: 

dfmM = 

a dt ( + \ ( 6 °) 

r — = 0- 

dt 

As mentioned above, under the steady-state conditions in principle there is no limitations 
in the use of the rate equations, except the steady-state condition (jSOjl itself. 



D. The case of large Zeeman splitting. 

In the case of large Zeeman splitting, that is, when (|39j) does not hold, the final rate 
equations become more complicated. However, the derivation procedure, of course, is still 
the same: we assume the already mentioned conditions, then simplify optical coherences 
p^T. and p^ 9i from (p5fi|) , (|37j) , and substitute them in (p?5|) , (|5H|) . Thus we arrive to the final 
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rate equations, which now become more complicated than ([44)1 . ()45|). The definitions of T p 
and AE p also become different than those in (}4Tj) - (|5Uj) . All the above analysis still holds. 
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